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Abstract 

We study localization in two- and three channel quasi-lD systems using 

multichain tight-binding Anderson models with nearest-neighbour interchain 

hopping. In the three chain case we discuss both the case of free- and that of 

periodic boundary conditions between the chains. The finite disordered wires 

are connected to ideal leads and the localization length is defined from the 

Landauer conductance in terms of the transmission coefficients matrix. The 

transmission- and reflection amplitudes in properly defined quantum channels 

are obtained from S'-matrices constructed from transfer matrices in Bloch 

wave bases for the various quasi-lD systems. Our exact analytic expressions 

for localization lengths for weak disorder reduce to the Thouless expression for 

ID systems in the limit of vanishing interchain hopping. For weak interchain 

hopping the localization length decreases with respect to the ID value in all 

three cases. In the three-channel cases it increases with interchain hopping 

over restricted domains of large hopping. 
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I. INTRODUCTION 



A wire is a topologically one- dimensional system whose width of the order of the square 
root of the cross-sectional area A is much smaller than its length L i.e. \fA « L. In a 
thin wire the motion of electrons in the transverse direction is quantized. The corresponding 
transverse eigenstates below the Fermi level define a finite number, N oc yf^ (with \p, the 
Fermi wavelength), of quantum channels for transmission of electrons across the wire. The 
starting point of the present work is the Landauer two-probe conductance formula |T||J] 



g = — Tr(it + ) , (1) 

which describes current transport in a disordered wire. Here t is the so-called transmission 
matrix of the iV-channel system: 



tn ti2 ■ ■ ■ tlN 

hi 

t = . (2) 

\tNl tm ■ ■ ■ tNN J 

An incoming wave from channel j of an ideal lead at one end of the disordered wire (of 
length L) has a coefficient \tij\ 2 for transmission into channel % of the lead at the other end. 

The correctness of the description of the conductance in terms of transmission channels 
has received striking experimental confirmation || in the special case of perfectly transmit- 
ting channels, where ([!]) reduces to 



9 = -N . (3) 

These studies relate to quantum point contacts in the form of narrow conducting two- 
dimensional strips whose width, hence the number of discrete transmitting channels, may 
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be varied by varying externally applied gate voltages. Conductance steps corresponding 
to increasing values of N in (^) are clearly observed || . A further interesting feature of 
these experiments is that they offer the possibility of realizing physically not only purely 
one-dimensional systems (N = 1) but also few-channel systems (N = 2, 3, . . . ) such as those 
studied below. We recall that recent discussions of many-channel mesoscopic systems have 
been essentially restricted to the case iV >> 1 |3| i.e. real (metallic) wires whose widths are 
much larger than the Fermi wavelength. 

The behaviour of the conductance in a disordered wire depends strongly on its length 
relative to the localization length L c . In an infinitely long disordered wire of a given cross- 
section all eigenstates are expected to be localized just like in the truly one-dimensional 
case. Thouless has indeed shown 0|| that for L > L c , where 

L C ~N£ , (4) 

I being the mean free path, the conductance should fall off exponentially as e~ L ^ Lc , which 
is a clear manifestation of localization. On the other hand, one the scales of lengths in the 
domain 

£ < L < L c = N£ , (5) 

the eigenstates appear as being delocalized. In fact, it was shown later |7| that for lengths 
in the range (5) there remains a number N e g ~ of independent ballistic channels leading 
to metallic Ohm's law behaviour for the conductance (|]). Such a diffusive quasi-metallic 
domain does not exist for truly one-dimensional systems where L c = L lc ~ £. Dorokhov || 
has developed a detailed scaling analysis || of localization in a multichannel wire, in which 
he calculates L c in terms of a phenomenological mean free path entering as input via Ohm's 
law at short scales. His final result, which is valid for weak disorder, is 
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L C = (N + 1)£ , (6) 

which coincides with (^) for iV >> 1 and suggests that L c is not directly proportional to 
the number of channels for few channel systems. In particular, for one-dimensional systems 
it yields 



L lc = 2£ . (7) 
We recall that (0) coincides with the result obtained long ago by Thouless Jl0| from kinetic 



transport theory. The same result follows also from the quantum series composition law for 
conductances g = \t\ 2 for ID conductors, assuming Ohm's law to be valid at short length 
scales [|llj . On the other hand, we note that for two-channel systems (|6|) yields L c = L 2c = 3£ 



which is comparable to the result L 2c = 2£(1 — obtained by Dorokhov [12j in a different 
scaling treatment. 

Shortly before the development of the scaling theory of localization for ID conductors [|9] , 
Thouless derived his well-known analytic expression for the quantum localization length in a 
tight-binding linear chain for weak disorder [13]. The purpose of the present work is to derive 



similar exact microscopic expressions for the localization length in quasi-one-dimensional few 
channel systems, specifically for two- and three chain systems for weak disorder. This study 
is of interest in several respects. For example, via the dependence of the localization length 
on the energy across the energy bands of the pure systems, it provides a first principles 
quantum proof of the fact that all states in the quasi-one-dimensional systems are localized. 
It is also relevant for experimental situations e.g. for discussing the quasi- metallic domain in 



few channel systems which are encountered in semiconducting microstructures ]14| or may 
be fabricated artificially, as in the quantum point contact system discussed above. 

In II. a we define the tight-binding two- and three chain systems for modelling few chan- 
nel wires. In II. b we describe our methodology for studying the localization length. It 



consists in constructing successively transfer- and transmission-reflection matrices for the 
tight-binding systems by generalizing well-known methods for ID systems. It further relies 
on general results concerning the existence and the properties of a Lyapunov exponent (in- 
verse localization length) describing the asymptotic exponential decay of the conductance 
(0) in mutichannel systems [T5|. In Sect. Ill we present the details of our calculations 
leading to the final analytic expressions for the transfer and scattering matrices for the two- 
and three channel wires for weak disorder. For the case N = 3, we obtain different results 
for open- and for periodic boundary conditions which correspond to packing the chains on 
a plane and on a cylindrical surface, respectively. The final analytic expressions for the av- 
eraged transmission- and reflection coefficients and for localization lengths are discussed in 
Sect. IV. For clarity's sake some details of these calculations are relegated to an appendix. 



II. FEW-CHANNEL WIRES AND LOCALIZATION 

A. Anderson models in channel bases 

We describe two (N = 2)- and three (N = 3)-channel wires by Anderson models for two- 
and three coupled chain systems, respectively. The two chain Anderson model consists of 
parallel linear chains of Nl disordered sites each (of spacing a = 1 and length L = N^a) 
connected at both ends to semi-infinite ideal (non-disordered) leads. It is defined by the 
tight-binding Schrodinger equation which we write in the matrix form 



{ i u i \ 



E — Sin —h 

■h E — E2n 



/ A 



(8) 



where the <p l m denote the wave-function amplitudes at sites m on the chain i, h is a constant 
matrix element for an electron to hop transversally between a site n on chain 1 and its 
nearest-neighbour site n on chain 2. The site energies Ei m are random variables associated 
with the sites 1 < m < Ni of the disordered chain i, and Ei m = on the semi-infinite ideal 
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chains denned by the sites m > Nl and m < 1, respectively. The above energies, including 
E, are measured in units of the constant hopping rate along the individual chains. 

The coupled three chain (N = 3) Anderson model is defined in a similar way by a set 
of tight-binding Schrodinger equations, whose actual form depends, however, on interchain 
boundary conditions. For free boundary conditions, which correspond to arranging the 
parallel equidistant chains on a plane, the tight-binding equations are 



\ 

















[pi) 



E - e ln -h 

-h E-e 2n -h & ■ (9) 

-h E- e 3n j 

with the sites in the disordered sections of length L = Nia and in the semi-infinite ideal 
chain sections labelled in the same way as in the two chain case. On the other hand, in 
the case of periodic boundary conditions which correspond to equidistant linear chains on a 
cylindrical surface the Schrodinger equation is 



/ i j_ i \ 

<Pn+l + <Pn-l 



+ <#u 



\P 3 n+l + J 



\ 



E — Sin —h —h 

-h E-e 2n -h ifl (10) 

—h —h E — e 3ri 

As discussed in Sect. I, a quasi-one-dimensional wire is described by a collection of 
independent channels for wave transmission. Microscopic models for two- and three channel 
wires are obtained from the systems of tight-binding equations (]8|-|iOD by diagonalizing the 
interchain coupling terms in the equations describing the ideal leads. This indeed leads to 
independent quantum channels for the leads defined by new amplitudes bases 



( '■■ ) 




( '■ \ 
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= u- 1 




V ; > 







in which the non-random parts of the matrices on the r.h.s. of (p|-|T0|) are diagonal. For the 
multichain system above we obtain, respectively. 
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-h E 
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-h E -h 
-h -h E 
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\ 



E — h 
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( 



U '=2 



\ 



1 y/2 1 ^ 
V2 -y/2 

1 -y/2 1 j 

( E - y/2h 
E 
, E + \/2hl 



\ 



1 1 
1 -1 -1 



E — 2h 
£ + 
E + h 



(8.a) 



(9.a) 



(lO.a) 



Note that £7o an d C/' are unitary, while U" is not. 

Finally, in the channel bases defined by (8a- 10a) the tight-binding equations (PlfTOQ read, 
respectively, 



€+1 + ^n-l 



E — h — \{Si n + £2n) ^(^2n — ^ln) 

^(^2n — £ln) E + h — \{Si n + £2r. 



(8.b) 
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^+1 + V£-1^ 



£- v/2/i- i(£l„ + 2£ 2n + £ 3r . 

"If ( £ 3n — ^ln) 
— li^ln — 2£?2n + ^3n) 



.E — 2/l — |(£i n + £ 2n + ^3r, 
— |(2£ln — 2e 2 n — ^3n) 
~ 2£2n + ^3n) 



"lf( £ 3n — £ln) 



~~ j( £ ln ~~ 2e 2 n + ^3r. 



V2 



( £ 3n — 



4 



(^3n £ ln 



(9.b) 



-Uei n ~ £3n) 

E + h - l(e ln + e 3 , 
i 



3 ( £ 2n — £ 3n, 



£ ln — £ 3ri) E + h — |(2e 2rl + £ 3n ) y 



(lO.b) 



which constitute our starting point for deriving transmission- and reflection matrices of the 
disordered wires in Sect. III. It is seen that the similarity transformation of the disorder 
matrices by the [/-matrices leads to interchannel coupling in the disordered sections, 1 < 
N < N L . 



B. Localization from conductance 

In Sect. IV we will calculate the localization length in the above multi-channel wire 
models from the rate of exponential decrease (Lyapunov exponent) of the conductance ([[]) 
for the large L [ . The transmission matrix in (Q) will be found by constructing a transfer 
matrix which transforms propagating waves in the multichannel leads (defined by (8.a-10.a) 
on the left side of the disordered wire into corresponding propagating waves on the right 
side. The transfer matrix for the wire of length L is expressed, as usual, as a product of Nl 



transfer matrices for small sections enclosing only the nth site of each one of the channels. 
The calculation of the localization length rests on theorems of Oseledec |TjJ and of Tutubalin 
and Vister |T7[ on the properties of products of a large number of random matrices. Indeed, 



employing these properties Johnston anf Kunz |T5j have shown that the Lyapunov exponent 
7 exists for the conductance (|I|) and is a self-averaging quantity referred to as the inverse 
localization length. It is defined by the relation 

where (...) denotes averaging over the disorder (i.e. the random site energies in the Anderson 
model) [Tj|. It follows that the asymptotic distribution of the conductance is log-normal. 

III. DETAILED ANALYSIS 

As indicated above the transmission matrices of the form (Q) for the quasi ID-disordered 
systems above will be obtained from transfer matrices for wavefunction amplitudes defined 
from equations (8.b-10.b), respectively. The construction of these transfer matrices proceeds 
in two steps. First, we will define transfer matrices for thin slices enclosing only one site n 
of each chain in a disordered wire, in a Bloch plane wave basis. Next the transfer matrix of 
a whole wire of length L = N^a will be obtained as a product of the transfer matrices for 
the Nl individual slices composing the wire. We will express it analytically to lowest order 
in the effect of a weak disorder. Finally we obtain the form of the tight-binding equations 
describing the transfer of Bloch wave amplitudes across a whole disordered wire of length 
L, which we then cast in the form of scattering equations in order to identify microscopic 
scattering matrices of the form (|2]), in terms of elements of the transfer matrices. We recall 
that the transfer matrix method is well-known in the study of one-dimensional disordered 
systems [PT5[— . Here we generalize and adapt it in the case of few channel quasi- ID systems. 
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A. Transfer matrices 



Transfer matrices X 0n , X' n and X n " for thin slices including a single site n per chain 
of the quasi-lD systems described by (8.b-10.b) are defined by rewrtiting these equations 
respectively in the forms 











— Xq u 
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where 
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< r' n 



E-ri -1 v' n 

10 

v' n E + V2h-/j' n -1 



\ 



2e 







V2, 











2n 







/ 



(15) 



(15.a) 
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1 1 1 

A*n" = -(^ln + ^2n + £3n) , ^ra" = g(2£l n + ^3n) , #n" = ~{^2n + ^3n) 

^n" = — ~(^ln _ ^3n) , ^n" = — g(^2n _ ^3n) , 

1 1 

«n" = — -(2£ln — £2n — ^3n) , /V = ~~ ^£2n + ^3n) 



(16) 



;i6.a) 



The study of disordered wires in terms of reflection and transmission properties of plane 
waves requires determination of plane wave bases in which the transfer matrices for slices n 
in the leads are diagonal. Such bases are provided by the Bloch wave solutions for the leads 
which are defined by 
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oo 



,±ifci„/,l 



(17) 
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(18) 



where Xoo, X' Q and Xo" denote the transfer matrices for the leads given by (|TJ]-[n^) with 
£in = ^2n = ^3n = 0, respectively. The wavenumbers, k h are defined in terms of the energy 
E by the eigenvalues of Eqs. (|17]) and (|T8|) , respectively By solving for the eigenvalues we 
get successively 

2 cos k\ = E — h , 

2 cos k 2 = E + h , (17.a) 

for the two-channel system, 

2 cos fci = E - V2h , 
2 cos tv2 = E , 

2 cos A; 3 = £ + y/2h , (18. a) 

for the three-channel system with free boundary conditions, whose leads are described by 
Xq, and, finally, 

2cosA;i = E-2h , 

2 cos k 2 = 2 cos k 3 = E + h , (18.b) 
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for the three-channel model with periodic lateral boundary conditions. The eigenfunctions 
in the leads at energy E obtained from (|T2| - |T3| ) and (|T7|-|T8|) are of the form 



3 ~ giinfcj 



n,± 



(19) 



where we choose the wavenumbers kj,j = 1,2,3 to be positive, < kj < tt, so that these 
functions correspond to plane waves travelling from left to right and from right to left, 
respectively. 

The transfer matrices for single site slices in Eqs (|12"1-|13|) for the leads (i.e. for Ei n = 0, 
n < 1 or n > N) are diagonalized in the bases of the Bloch plane wave states (|T9"D. The 



diagonalization matrices which are formed by the eigenvectors of (U7- 18), are of the form 



In 



A; 



^ikj g ikj 



y/ 2i sin kj 



1 J 



for the two-channel quasi- ID model, and 



(20) 



V 



O A 2 O 
O A 3 



(21) 



with Aj defined as in (|2"0"|), for the three-channel models. The wavenumbers kj defined by 
(17a) for the N = 2 case and by (18. a) and (18b) for the N = 3 case with free- and periodic 
boundary conditions, respectively. After finding the inverses of Vq and V and performing 
the simularity transformations of X 0n , X' n and X n " by Vq and V, respectively, we obtain 
the desired transfer matrices in the Bloch wave representation of the disordered wires. In 
the two-channel case we find 



12 



Xo n = V 1 X 0n V 



where 



e ikl (l +ia ln ) ie~ ikl a ln 
-ie ikl a ln e- ikl (l-ia ln ) 
-i e ikl b r , 



i e ik2 b r , 



i e 



-i e ikl b n e ik2 (l + ia 2n ) ' ' ""-a 2 



-i e ik2 b ri 
i e-' lk2 b n 

—ik2. 



i e %kl b r 



i e tkl b r , 



■ie ik2 b n e- ik2 (l-ia 2n )J 



din 



Sir. 



&2n Sin + £ 2r . 
, <>2n 



4 sin ki 

£2n — £ln 



4 sin k 2 



4^ sin ki sin k 2 ' 



(22.a) 



(22) 



are real quantities and hi, k 2 are defined by (17a). For the three-channel systems we write 
the final transfer matrices X' n and X" n , in terms of a generic matrix 



Z n — 
( 



e ikl (l + ia ln ) 


i e lkl a ln 


i e lk2 c n 


i e~ ik2 c„ 


i e lk3 g n 


1 e-^g n 


-i e lkl a ln 




-i e lk2 c n 


-i e~ ik2 c n 


-1 e ik *g n 


-1 e-^g n 


ie ik ^f n 


ie- ik ^f n 


e ik *(l + ib 2n ) 


i e- lk *b 2n 


i e lk3 d n 


i e- ik3 d n 


-% e^fn 


-i e"^/„ 


-i e^b 2n 


e- ik2 (l-ib 2n ) 


-i e tk3 d n 


-i e~ ik3 d n 


i e lkl p n 


i er %kl Vn 


1 e lk2 q n 


1 e~ ik2 q n 


e ik3 (l + ia 3n ) 


i e- lk3 a 3n 


-i e lkl p n 


i e~ tkl p n 


-i e^q n 


-i e- lk2 q n 


-i e lk3 a 3n 


e- ik *(l-ia 3 



, (23) 



V 



where ai n , a 3n , b 2n , c n , g n , f n , d n ,p n and q n are real quantities. In the case of free boundary 
conditions we obtain X' n = V^X^V = Z n with 

_ Eln + 2^2n + ^3ra _ £ln + ^2n + ^3n 

O-ln — — ; , «3n — " 



b 2n = 



8 sin k\ 

gin + e 3n 
4 sin k 2 



fn 



8 sin fc 3 

V2(£ln ~ ^3n) 

8Vsin k\ sin k 2 



, _ V / 2(£ln — ^3n) _ £ln — ^2n + ^3n 

7i / ■ , ■ , j 9n Pn 



(24) 



8vsi n k 2 sin A; 3 ' 8\/sin k\ sin /c 3 

Here /c, ki, /c 3 are defined by (18. a). On the other hand, for periodic boundary conditions we 
find X"„ = V^X^V = Z, n where 
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E\ n + E2n + £ 3n 2e 2 n + £ 3n 

a ln — n ) 



6 sin A;! 6sinA; ; 



2 



, _ 2£i n + E^n Sin ~ £ 3n £ 2n ~ £ 3n 

"2ra « ■ i ' ^ n 7> / ■ 7 t : ' 9n 



6sinfc 2 6A/sin fci sin k 2 6\/sin k[ sin /c 2 

, _ ^3n — ^2n n _ ^ £ ln ~ £ 2n ~ £ 3n 
8(1 — T, : ~j j In 



6sin/c 2 6Vsin k~i sink 



2 



~ £ ln + 2£2n — £ 3n £ 3n ~ £ ln , or N 

6\/sm ki sin fc 2 6 sin fc 2 

where &i and k 2 = k 3 are now given by (18b). Note, in particular, the diagonalization of the 
transfer matrices for the leads in the plane wave bases shown in 

Finally, we determine the transfer matrices for the disordered wires of length L = N^a 
in terms of the transfer matrices of the individual thin slices n. As shown by iteration of the 



transfer equations (|T7| - |T8|) rewritten in the Bloch wave basis above, the matrix transferring 
an incoming wave at site n = just outside a disordered wire to the site Nl + 1 just beyond 
its other end is given by a product of transfer matrices of the form 

Nz, 

Y L = l[Y n , (26) 

n=l 

where Y n and Y L stand for the three pairs of transfer matrices X 0n , X 0L , X' n , X' L and 



respectively, which are associated with the wire models above ||19fl . 
We shall evaluate the transfer matrices of the disordered wires for weak disorder to 
linear order in the random site energies. On the other hand, for our explicit calculations of 
averages over the disorder below, we assume the site energies to be independent gaussian 
random variables with zero mean values, and correlation 



( £ in £ jm) £ Q^i,j^n 



(27) 



In this case the site energies corresponding to different slices n in (|2q ) are uncorrelated so 
that it is indeed sufficient to restrict the expansion of the latter expressions to first order 
in the site energies for determining averages to order e§. The transfer matrices Y n given by 



23) are sums, 
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(28) 



of a zeroth order diagonal matrix y(°) independent of the site energies (transfer matrix of 
the leads) and a matrix which is linear in the energies Sj n . By inserting ( p8[) into (|26|) we 
obtain to first order 



Y L = (Y^) Nl + ^(y(o))™-iyW(yW)^-m + . 



(29) 



m=l 



Next we insert the slice matrices Y^ and Ym^ from (p2| - ^3D for our various quasi-lD systems 
and obtain successively: 

X 0L = diag(e ifclA \ e- iklNL , e ik2N \ e ~ ik2NL ) 

^ ia lm e lkl ui ia lm e~ lkl sl -ib m e lk2 w 2 i -ib lm e~ lk2 vl^ 

-iai m e tkl si -ia lm e~ tkl u\ ib m e lk2 v 21 ib lm e~ %k2 w\\ 

-ib m e lkl w 12 -ib m e~ tkl v* 12 ia 2m e lk2 u 2 ia 2m e~ tk2 s 2 

y ib m e tkl v 12 ib m e~ lkl w* l2 -ia 2m e lk2 s 2 -ia 2m e~ tk2 u* 2 J 



Nr 



m=l 



for the two-channel case. Here 



Sj = e 



i(N L -2m+l)kj 

, Uj 



i(N L -l)k 3 



Vij = e 



i(N L -m)ki-i(m-l)kj w .. — e 



i(Ni,—m)ki+i(m—l)kj 



N L 



where k± and k 2 are defined by (17a) 

X' L = diag (e'W^-'W^'W^-W^'W^-Wi) + ^ 

m— 1 



iai m e lkl ui 


ia lm e ikl sl 


ic m e tk2 w 2 i 


ic m & ik2 V2i 


ig m e lka w 31 




-iai m e lkl si 


-ia lm e~ ikl u\ 


-ic m e lk2 v 2 i 


—ic m e ik2 W2i 


-ig m e lk ' 3 v 31 


-ig m e- lk3 w* 31 


if m e lkl w 12 




ib 2m e lk2 u 2 


ib 2 rae~ %k2 S*2 


id m e lk ' 3 w 32 


,7/ p -ik3„,* 
ta m e U32 


-if m e lkl v 12 


*/m e ikl Wi2 


-ib 2m e lk2 s 2 


-ib 2m e- ik2 u* 2 


-id m e lk3 v 32 




ip m e lkl wiz 


iPm e iklv 13 


iq m e tk2 w 2 3 




ia 3m e ik3 u* 3 


ia 3m e~ ik3 s* 3 


-ip m e lkl v 13 


-ip m e- ikl wl s 


-iq m e' lk2 W23 


-iq m e~ ik2 w* 23 


-ia 3m e ik3 s 3 


-ia 3m e~ ik3 u* 3 



(30) 



(31) 



(32) 
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for the three-channel case with free boundary conditions where k\,k2,ks in the definitions 
(p4]) and fl3~l"|) are given by (18. a); the transfer matrix for the three-channel model 
with periodic boundary conditions is given by (|32|), using the definitions (|25| ) of the slice 
parameters ai m , a 3m , & 2m , c m and the definition (18b) of the wavenumbers ki, &2, k 3 . 



B. Scattering matrices 

The scattering of plane waves (reflection and transmission) at and between the two ends 
of the random quasi- ID systems is governed by the S'-matrix, 



S 



r + t 



t ++ f, 



(33) 



where 



and 
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/ r ±T 
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r ±T 
'21 
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(34) 



(35) 



: / 



Here ^{t^ ) and r 



denote the transmitted and reflected amplitudes in channel i 



when there is a unit flux incident from the left (right) in channel j. Left to right- and right 
to left directions are denoted by + and -, respectively. The S'-matrix expresses outgoing 
wave amplitudes in terms of ingoing ones on either side of the quasi-lD disordered wire via 



the scattering relations [11,21 
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Here I and I' (0 and 0') denote ingoing (outgoing) amplitudes at the left and right sides of 
the disordered region, respectively. It follows from current conservation that e.g. for a unit 
flux which is incident from the right in channel i one has 

N 

Ed *7rr + k~+ 1 2 ) = i . (37) 

Likewise, one has also 

E(l t%+ ? + I rt- | 2 ) = 1 . (37.a) 
i=i 

Our task is now to derive miscroscopic realizations of the ^-matrix in terms of the transfer 
matrices (j3~C|) and (^) describing transfer of Bloch waves across finite quasi- ID disordered 
systems. Let us remark that this may not always be possible as is seen here in the case 
of the three channel periodic model. Indeed, in this case we are able to identify proper 
transmission and reflection amplitudes obeying the symmetry relations fl3?D and (37a) only 
when assuming the random site energies on chains 1 and 2 to be identical, 



£\n = e 2n , n = 1,2,... N L , (38) 

rather than allowing the energies at all pairs of sites of the quasi- ID system to be uncorre- 
lated, as in (|27|) . In the absence of the correlation ( j38|) the obtained transport amplitudes 
cannot be identified as actual reflection and transmission amplitudes of Bloch waves. 

If the amplitude at the nth site in a channel j corresponds to a Bloch wave ipi^ = e mfcj 
(equation (|I9D) then the j,n and j, n — 1 components of wave amplitude vectors, 
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( '■■ ^ 

w- 1 Wn 

(with W = V or V), being transferred by the n-th slice have values a e mfej and 0, respec- 
tively; on the other hand if ip 3 n = e~ mk ^ then the j, n and j, n — 1 components of the above 
vectors are and a'e~ mkj . Hence, in accordance with our notation above for transmission 
and reflection matrix elements we denote the amplitudes at site j, n and j,n — l respectively 
as o^ n „i and Q>J n -i since they correspond to amplitudes being transferred by the n-th slice 
and propagating in the kj and —kj directions, respectively. The transformed amplitude 
vectors 



and 



V : J 



w 1 



( : \ 

V ; J 



are thus rewritten, respectively, as 
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1 11 \ 


/ + \ 

fl l,n-l 








( + \ 

a l,n 




a l,n-l 








a l,n 


r n = 


= a 2,n-l 


and W^ 1 






a 2,n 


^-1 


a 2,n-l 








a 2,n 



(39) 



V : / V : / V : / V : / 

Using a similar notation for wave amplitudes transferred from n = to n = Nl across a 
disordered wire of length L = N^a, the wave transfer equations in the Bloch representation, 
obtained by iterating (IT^-O) read 



( <L ) 




"l,0 


<h,L 


— x 0L 


a l,0 


a 2,L 




a 2,0 






\ a 2,0 ) 



(40) 



and 



w 




u l,0 


a l,L 




°1,0 


a 2,L 




a 2,0 


a 2,L 




a 2,0 






a i,o 


\ a 3,L j 




\ a 3,0 J 



(41) 



where X 0L is given in (|30D and Y L stands for X' L and X" L in ( |32|) with parameters defined 



by (18. a) and (|24f) and by (18. b) and (|25f) , respectively. 

In order to derive the S-matrices for our two- and three channel wire models we first 



rewrite (JD-JT) in the forms of equations involving outgoing amplitudes on the left side and 
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incoming ones as the right side, as in fl36|). In the notation of (f40"D and (|4lD we have e.g. for 
N = 3 



b i,o 



'2,0 



(/') 



1,L 



A 2,L 



(O) 



1,0 



l 2,0 



(O0 = 



1,L 



2,L 



(42) 



\ a 3,0 / 



\ a 3,i / 



\ a 3,0 / 



\ a 3,L J 



and so we rearrange ( |40"| - f41"| ) in the form 









(A 


^2,3 














W 



(43) 



where the pairs of matrices A 2 and B 2 and A 3 and B 3 correspond to two- and three channel 
cases, respectively. Using the notation (X 0L )ij = Xij and (Y L )ij = Y^- for the matrix 
elements in ( |3T)D and respectively, we find 



( 



A, 



X%2 ~ X\± 1 0^ 
—X22 —X24 
— X 32 — X34 1 

-x 42 -X44 oy 
\ 



s 2 



\ 



-Y 12 


-Y u 


-Yl6 


1 








—Y22 


—Y 2 4 


—Y 2 e 











— Y32 


-Y M 


— ^36 





1 





—Y42 


— 144 


— ^46 











-Y 52 


— Y54 


— ^56 








1 


— Yq2 


— Yq4 


— ^66 












X u X 13 
X21 X23 —1 

-^31 X 33 

yX 4 i X 43 -1J 
I 



(44) 



ni ^13 n 5 

y 21 y 23 y 25 -1 

Y 31 Y 33 F35 

F 4 i Y 43 F 45 -1 

*si n 3 155 

>61 *63 *65 0-1/ 



(45) 



The S'-matrices for the two- and three channel cases are those given by S = A2 B 2 and 
S = A 3 ~ 1 B 3 , respectively. After inverting A2 and A 3 we finally obtain, successively for the 
two- and three channel cases: 
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/ 



where 



Si 82 X44 — X24 

S 3 S4 —X42 X 2 2 

X U S + X 2 i5 5 X 13 S + X 23 5 5 —S 5 -S 6 

X 31 S + X 2 i5 7 X 33 S + X 23 5 7 -S 7 -5 8 

+X4\5s +X 43 S 8 



\ 



S = X22X44 — X24X42 , S\ = X24X41 — X44X21 , S2 = X24X43 — X44X23 , 

S 3 = X 42 X 2 l — X 41 X 2 2 , O4 = -^42^23 — -<^22^43 j ^5 = -^42-^14 _ X U X 4 4 , 

Sq = X12X24 — X22X14 , S 7 = X42X 3 4 — X32X44 , Ss = X 3 2X24 — X22X 3 4 , 

are second order subdeterminants of X 0L ; 



(46) 



(47) 



^ Si S 3 

yS*2 5*4 J 



where 



Si 
( 



1 

~ A 

-PiY 2 i - (3 4 Y 4 i - [3 7 Y 6 i -(3iY 23 - ^ - f3 7 Y 63 -(3iY 25 - (3 4 Y 45 - (3 7 Y 65 

P2Y21 + fcY 4 i + foYei P2Y23 + fcY 43 + /3 8 F 63 (3 2 Y 2b + fcY 45 + /5 8 F 65 

~/^3^21 — PqY 4 i — P9Y51 —f3 3 Y 23 — (3qY 43 — (3gY 63 —f3 3 Y 2 5 — /?6^45 — ^9^65 



(48) 



(48.a) 



^Y n A + Y 21 A 1 + Y 41 A 4 + Y 61 A 7 Y 13 A + F 23 Ai + Y 43 A 4 + Y 63 A 7 Y 15 A + Y 2b A 1 + Y 45 A 4 + Y 65 A 7 ^ 



Y 31 A - Y ?A A 2 - Y 41 A R - Y M A s Y 33 A - Y 23 A 2 - F 43 A B - Y fi3 A 



'-21^2 - *41^5 



'61^8 



:s:s- 



[ 23*-V2 — ^43^5 



'63^8 



Y35A — Y" 2 5A 2 — I45A5 — Ye5 ^8 

F 5 iA + F 2 iA 3 + Y 41 A 6 + Y ei A 9 Y 53 A + Y 23 A 3 + F 43 A 6 + F 63 A 9 Y 55 A + Y 25 A 3 + Y 45 A 6 + Y 65 A g j 

(48.b) 



21 



;/ A ft fa ^ 



1 

A 



V 



-ft -ft -A 

ft ft ft 



(48.c) 



-Ax -A 4 -A 7 ^ 



5, 



1 

A 



A 2 A 5 A, 
-A 3 -A 6 -A 9 y 

which involve subdeterminants of second order of Yl, 



V 



(48.d) 



ft = ^46^64 — Y 44 Yq6 , ft = Y^Xm — Y^Yqq , ft = Y^Y 44 — Y42YQ4 , 
ft = ^24^66 — ^64^26 j ft = ^22^66 ~~ ^62^26 j ft = ^22^64 — ^62^24 > 

ft — Y44Y2Q — Y24Y4Q , ft = Y42Y2Q — Y22Y4Q , ft = Y24Y42 ~ Y22Y44 , (48. e) 

as well as third order subderterminants of A 3 which result from the minors of various ele- 
ments, 

Ai = min(A>) 24 , A 2 = min(A>) 25 , 

A 3 = min(A 3 ) 2 6 , A 4 = min(A 3 ) 44 , A 5 = min(A 3 ) 4 5 , 

A 6 = min(A 3 ) 46 , A 7 = min(A 3 ) 64 , A 8 = min(A 3 ) 65 , 

A 9 = min(A 3 ) 66 , and A = det A 3 . (48.f) 

The partial matrices Sj, j = 1, 2, 3, 4 in (f48"|) clearly correspond to reflection and transmission 
matrices in fl33|). 

The matrices ( f4~6|) and (^8],^8]a-ci) are the S'-matrix expressions for the two- and three 
channel quasi- ID systems in terms of characteristic quantum channel wavenumbers and of 
the tight-binding quantities in the transfer matrices (|30| ) and ( |32"D for weak disorder. Iden- 
tification of these expressions with (|33| - |3"5|) yields the transmission- anf reflection matrices 
for these quasi-lD disordered models under the proviso that the symmetry relation (P7|-37a) 
are obeyed. 
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IV. RESULTS AND DISCUSSION 



The results for the transmission- and reflection matrices of two- and three channel tight- 
binding wires are applied in this section for finding the averaged transmission and reflection 
coefficients associated with the various channels, for weak disorder. These results allow us, 
in particular, to explicitely check the symmetry property (|37|) in the two-channel case, as 
well as in the three-channel case with free boundary conditions. On the other hand in the 
three-channel case with periodic boundary conditions we show that ([37]) is obeyed if one 
restricts the disorder to a correlated site energy disorder with identical site energies Q38|) on 
chains 1 and 2 and independent random energies on chain 3. 

The results for the averaged transmission coefficients are used for obtaining the length of 
exponential localization from ( |TT| ) and ([I]). Our exact quantum expressions of localization 
lengths for weak disorder reduce to the well-known Thouless expression for a ID chain in 
the limit of vanishing interchain coupling {h — > 0). It is useful, before presenting our results, 
to briefly recall the derivation of Thouless' result from the transfer matrix approach Jl9 
For a single disordered chain of length L the localization length is given by 



i = - Jim (2iV i )- 1 (ln|t| 2 ) , (49) 

where | t \ 2 =\ t | 2 =| t ++ | 2 is the transmission coefficient which is related to the two- 
dimensional transfer matrices X n (in the Bloch wave representation) for thin sections en- 
closing the nth site by W~ 



\n=l / 22 

This relation follows by transforming the S'-matrix for scattering states into a transfer matrix 
whose 22-element is l/t . Expanding the transfer matrix for the whole chain, n^=i-^rn to 
first order in the uncorrelated site energies and performing the average in (fj5|), using (P7|), 
yields 
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which is Thouless' expression for the localization length in the tight-binding band, for weak 
gaussian disorder jl3| . 



For convenience of the following discussion for two- and three channel systems, the 
explicit forms of the transmission and reflection coefficients obtained by identifying the S- 
matrix (|33"1-|35|) successively with flUf ) and ([48], 48.a-d) and replacing the transfer matrix 
elements entering in these expressions by their explicit forms in ( |30| - |3"2| ) are given in the 
appendix. 

A. Two-channel wires 

By averaging the partial transmission- and reflection coefficients given by (A.1-A.4) over 
the disorder, using (22. a) and fl2"T|), we obtain successively 



(| (lT |. > = 1 _^( * ■ (52) 

8 V sm h smkismk 2 J 



(l I') " 1 " ^ (— Tg + sirifcl 2 S mA; 2 J ' (53) 



(l^l 2 ) = ^Ar > ™ 

8 sin k\ 



(k 2 - 2 + i 2 } = £4r > <™> 

8 sin k 2 
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12 



2 ) = (I ^21 + I 2 ); 



NlsI 



(57) 



8 sin k\ sin k 2 

The expressions ( |52] - |57| ) are consistent with the symmetry relations (|37| ) resulting from 
current conservation. 

The inverse localization length for weak disorder is obtained by expanding (|ll|) to lowest 
order in the random site energies using (fj) and ([52|-[54]). It is given by 



L c L 0c 32 ysinfex sinfc 2 

for energies E restricted to the Bloch bands E = h + 2 cos 
exact expression for weak disorder reveals three important 

1. It proves miscoscopically that all states in the Bloch energy bands of two-channel 
quasi- ID disordered systems are localized. 

2. In the absence of interchain hopping (h = 0) it reduces to the localization length ( |5"T| ) 
for a ID chain described by the Anderson model. For weak interchain hopping (|58"D 
becomes 




hi, and E = —h + 2 cos k 2 . This 
properties: 



8 sin k 



1 + 



h 2 



4 sin k 



[1 + 3 cot 2 k) +0(h 4 ) 



(59) 



which is valid for energies sufficiently close to the band centre of the h = energy 
band, E = 2 cos k. This shows that a weak interchain hopping enhances localization 
in comparison to the purely ID case, i.e. Lq c < £. 

3. For large interchain hopping rates, i.e. | h \»\ E \ (where E is of the order of the 



fermi energy) we have sinA;i i2 ~ y/l — h 2 /A , | h | /2 < 1, which yields 
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where £o = ^ s the ID localization length ([51]) at the band centre. Thus at 

large hopping rates (with | h ]< 2) localization is also enhanced in comparison to ID 
localization. 

B. Three-channel wires 

1. Free boundary conditions 

The evaluation of the disorder averages of the transmission- and reflection co- 
efficients in (A.5-A.6), using (|27|) and the explicit expressions in (A.7-A.9) with the 
tight-binding parameters yields, to order e^, 



„ x N L e 2 n ( 3 6 4 

( t — ) = 1 ( I I , (61) 

11 32 ysin 2 /?! sin &i sin £3 sinfcisinfc 2 



/| t — ,2\ _ 1 N L ej ( 4 | 4 | 4 , 
22 32 V sin 2 k 2 sin fc 2 sin k 3 sin /ci sin fc 2 



/. ,- ,2v _ ! _ ^£o / 3 6 4 . 

Vl % I / — 1 32 I sin2 h + gin ^ gin ^ + gin k2 gin fca ) > 



(I t" I 2 ) = (I t~ I 2 ) = 2 ^ Lg ° - , (64) 

Vl 12 1 7 \! 21 1 / 32 S in fci sin k 2 ' v 7 
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(I T.+ R = , (67) 



32 sin k\ 



(\^ + \ 2 ) = 4^r > ( 68 ) 

32 sin re 2 



(IrJIIHkJf) ^; , (71) 

sm ki sin /c 3 



(I r 2 - 3 + | 2 ) = (| r 3 -+ | 2 ) . 2 f Leg , • (72) 

sm K2 sm K3 



Again, the expressions (^TJ-[T^) obey the current conservation property (|3"7 ) for the 
three channel case: the reduction of the intra-channel transmission coefficients due to 
scattering by the disordered wire is exactly compensated by the occurence of inter- 
channel transmissions and by reflections. 

For the inverse localization length we obtain, from (0),(ffl) and (jBTj-BB]), 
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1 



+ 



+ 



64 ^sin 2 hi 3 sin 2 k 2 sin 2 k 3 



+ 



+ 



+ 



3 sin k\ sin k 2 3 sin k 2 sin k 3 sin k\ sin 



(73) 



where ki,k 2 , k 3 are denned by (18. a). Like (|58|), this expression is exact to order for 
weak disorder. It demonstrates that the eigenstates in the Bloch bands of the three 
channel quasi-lD system with free boundary conditions are localized. It too reduces 
to the Thouless result ([51]) in ID for vanishing interchain hopping. For small values 
of | h |, in particular for | h \« 1 for E — > 0, (|73|) reduces to 



1 



L c 8 sin k 2 



h 2 



sin k 2 



1 + i cos 2 k 2 + i cot 2 k 2 ) O('ir) 



and for large values (h 2 » E 2 /2) restricted to | h \ < 2 it becomes 



1 



8 2-h 2 



(74) 



(75) 



which shows, in particular, that L c is increased with respect to the ID value at the 
band centre (£o) for | h \ (in units of the interchain hopping rate) less than 1. Thus, the 
domain of interchain hopping rates in which ([73]) leads to localization lengths larger 
than £ is defined by 



E |« V2\h\<V2 



(76) 



2. Periodic boundary conditions 

We first note that in this case the S"-matrix constructed from the transfer matrix 
Xl" in Sect. Ill is unitary (i.e. ([37]) and (37a) are verified) only in the case where ( |38D 
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is obeyed i.e. when chains 1 and 2 have identical random site-energies while the ran- 
dom site energies on chain 3 are arbitrary. This may be readily seen by expressing the 
averages of the transmission- and reflection coefficients in (A. 5) and (A. 6) to second 
order in the site energies in terms of the averages of (A. 7- A. 9) assuming that site ener- 
gies belonging to sites m^non the same chain or on different ones are uncorrelated 
i.e. (Ei m Ejn) = 0. In this way we find 

E ((I t-f | 2 ) + (| r 4 r+ | 2 » = 3 + 2N L [«c 2 > + (ft) - 2(c n f n )) 

+(K> + (ql) - 2(d n q n )) 



+ ((9 2 n) + (pl)-2(9nPn))] • (77) 

Now, for this relation to be compatible with current conservation we require 



c„ = f n , d n = q n and g n = p n , (78) 

and from the definition (p5|) it follows that the equalities ( [78|) are fulfilled with (|38D 
i.e. for identical site energies in the chains 1 and 2. Only under this condition does 
the S"-matrix for the periodic 3-channel system represent a true scattering matrix 

The explicit second order expressions for the averaged transmission- and reflection 
coefficients obtained from (A.5-A.9) and (25), for random gaussian site energies ( |2"TD 



in the presence of the correlation (B3) are 



hl |2>_1 3V (sin 2 *! + sin/t lS m/fc 2 I ' (?9) 
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(I l 2 > = <l *5T l 2 > = <l *JT l 2 > = <l *5T l a > - 2 * ViL "" 1 



36 sin ki sin k 2 



(81) 



(l I > = (l & n 



2 . 2N L el 1 



36 sin 2 k 2 



(82) 



(ku + I 2 ) - BNl£2 ° 1 



36 sin fci 



(83) 



(I r 22 + | 2 ) = (| r 3 - 3 + | 2 ) 



(I rr 2 + T) 
(k 2 - 3 + r) 



,2-. /, — ,2-. 5N l4 1 
33 I 

<I'-Jf | 2 > = W | 2 ) = (I r 31 + f) - — ' 



36 sin 2 k 2 ' 



36 sin k\ sin fc 2 



(k 3 - 2 + l 2 ) 



2iV L£ 2 1 
36 sin 2 ko 



(84) 



(85) 



The inverse localization length associated with the conductance ([!]) obtained from 



91 32 is 



1 e 



14 



+ 



L c 216 V sin kx sin k 2 sinfcisinfc : 



(86) 



where k\ and k 2 are defined by (18b). This expression which proves localization in 
periodic 3-channel systems, reduces again to the ID result (|5"1"D for h = 0. For small /i 
it is (with E = 2 cos fc) 



1 



L c 8 sin 2 k 



h L 70 2 , 

1 H =— H cot 2 

2 sin 2 k \ 27 



(87) 
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which shows enhanced localization for weak interchain hopping. On the other hand 
for I h \»\ E I (EE) becomes 



5e 





L r 216 1 - h 2 



2 1 

(88) 



which leads to increased localization lengths, L c > £o, for h 2 < 22/27. Thus the 
domain where the localization length in the periodic 3-channel system is larger than 
the ID-value is defined by 



| E |«| h |< 0.906... . (89) 

Our exact microscopic results for two and three-channel systems indicate that 
in all case while in the three channel weak interchain hopping decreases the 

localization length from its ID- value. For strong interchain hopping a similar decrease 
persists in the two-channel case while in the three-channel cases the localization length 
increases with respect to the ID result over restricted domains of hopping parameters 



defined by (|75|) and (|jj). This suggests the possible existence of quasi-metallic do- 
mains for the three-channel systems in some ranges of length scales lying between 
their localization lengths and the ID localization length (of the order of the mean free 
path £ |TI| ) , in the considered ranges of large hopping. In contrast , we recall that for 
many-channel random wires (N » 1) the metallic domain exists over a wide range 
of mesoscopic lengths defined by (|5|) . 

V. APPENDIX 

For briefness' sake we only discuss explicit expressions for transmission amplitudes tjf 
and reflection amplitudes r~- + relating to outgoing waves to the left of the disordered sample 
in the scattering equations ( |3~6"D (see (|33D and (|42|)). 
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1. Two-channel wires 



Identification of ([33]) and fl4"6|) yields 



I v 12 I v |2 
, 2 I A 44 | | ,2_ I A 24 | 



y ll I I ^ 12 ' I "12 I 1 ^ | 2 ' 

I V 12 I v 12 

+ |2 A 42 , I , |2 A 22 /a -. \ 

*21 I = I g 12 aIld I ^22 I = I g p , (A- 1 ) 



and, to lowest order in the random site energies, using (|47|) and ([30 



|2 y 12 I -+ |2_| v |2 

A 21 I ) I r l2 I — I A 23 



11 I — I ^21 I 5 I 7 12 I — I ^-23 I ? 
21 + I 2 = I -^41 | 2 , I r 22 + | 2= l ^43 | 2 • ( A -2) 



^From the definition of 5 in fl4"?D and the explicit form of the elements X t j in ( |30| - |3~1~] ) , 
we obtain 

N L 

I 5 | 2 = 1 + 2^ [oimQin + a 2m a 2 n + 26 m 6 n cos(m - n)(fci - k 2 )\ , 

m,n=l 
1*22 | 2 =1 + E 

Qlm"'ln 5 I X44 | 2 = 1 
m,n m,n 

I -X12 I 2 = ^ai m a 2n cos2(m - n)ki , | X 24 | 2 = ^ & m &n cos(m - n){k\ - k 2 ) , (A.3) 



-+ 
11 


|2 


= J \m a \n cos 2(m — n)ki 

m,n 


-+ 

12 


|2 


— °m°n cos(m — n) (k\ + k 2 






m,n 


- + 
21 


2 


= & m &n cos(m - n) (A;i + A; 2 






m,n 


- + 
22 


2 


= y] «2mfl2n cos 2(m - n)fc 2 



(A.4) 



2. Three channel wires 

The identification of (48. c) with the transmission matrix i in ([34]) yields 
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I*u~ 


,2 IA 
1 A 


2 

12 ' 


^12 


|2 1 A 
1 A 


2 

2 ' 


^13 


,2 \fr 
1 A 


2 

2 ' 




^21 


,2 1 ^ 

" I A 


2 

12 ' 


^22 


1 A 


12 
2 ' 


^23 


,2 Ift 

1 A 


2 

2 ' 




^31 


,2 1 A 

1 A 


2 

12 ' 


^32 


,2 1 & 
1 A 


12 
2 ' 


^33 


,2 1 A) 
1 A 


2 
2 


(A.5) 



Similarly, the reflection amplitudes r~ + are obtained by identifying Si, in (48. a) with 
(PS]). For analyzing the reflection coefficients for weak disorder it suffices to find the 
amplitudes to linear order in the site energies. Such linear contributions are associated 
with f3 j-terms having a zeroth order contribution (of modulus one). From (48. a) and 
(48. e) we thus obtain 



1 r ll 


2 _ 


Y 2l 


|2 

i 


1 r l2 + 


2_ 


^23 


|2 

> 


1 r 13 + 


2_ 


^25 


|2 




1 r 21 


2 _ 


Yax 


|2 

i 


r 22 + 


2_ 


^43 


|2 


r 23 + 


2_ 


Y45 


|2 

i 




1 r 31 + 


2 _ 


Y 61 


|2 

) 


r 32 + 


2_ 


^63 


|2 


1 r 33 + 


2_ 


Y&5 


2 


(A.6) 



The explicit expressions of the | (3j | 2 and | A | 2 in (A.5) (defined in (48.e,f) and of the 
reflection coefficients (A.6) in terms of the transfer matrix elements in (|32]) are given 
by 

f3i\ 2 = 1 + ^[a 3m a 3n + b 2m b 2n + 2d m q n cos(m - n)(k 3 - k 2 )} , 

m,n 

/5 5 | 2 = 1 + ^2[a lm a ln + a 3m a 3n + 2g m p n cos(m - n)(k 3 - ki)] , 

m,n 

P 9 \ 2 = 1 + ^J[ai m a ln + b 2m b 2n + 2c m / n cos(m - n)(k 2 - ki)] , 

m,n 

p 2 \ 2 = ^ f m fn cos(w - n) (h - k 2 ) , | /3 4 | 2 = ^ c m c n cos(m - n) (k x - k 2 ) , 

(3 3 | 2 = ^ PmPn cos(m - n) (ki - k 3 ) , | (3 7 | 2 = g m g n cos(m - ra)(fci - k 3 ) , 

/3 6 | 2 = ^2q m q n cos(m - n)(k 2 - k 3 ) , | /3 8 | 2 = J^cf m d n cos(m - n)(k 2 - k 3 ) , (A. 7) 



m,n 



= 1 + E 

[almoin + a-zmO-zn + b 2m b 2n + 2g m p n cos(m — n) (fc 3 — /C4) 
+ 2<i m g n cos(m - n) [k 3 - k 2 ) + 2c m f n cos(m - n) (k 2 — k{)] , (A. 8) 
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r 



11 



r 33 



+ |2 



^ QlmQln , | r 22 + | 2 = ^6 2 m&2nCOs2(m - n)/c 2 , 

^a 3m a 3 „ cos 2 (m - n)fc 3 , | rf 2 + | 2 = c m c„ cos(m - n)(fa + k 2 ) 



m,n m,n 

-+ |2 



r 13 



„-+ 12 
23 



r 32 



2J SW0Vi cos(m - n)(ki + fa) , I r 21 + | 2 = 2J fmfn cos(m - n)(A;i + fa) , 

m,n m,n 

^ dnA cos(m - n)(fc 2 + fa) , \ \ 2 = ^p m p n cos(m - n){fa + fa) , 

m,n m,n 

^q m q n cos{m-n)(fa + fa) , (A.9) 



m,n 



where the double summations run from m = 1 to m — Nl and n = 1 to n = N^. 
The site-dependent tight-binding parameters dy, a 3 j, b 2 j, Cj, dj, fj, gj,Pj, qj in the above 



expressions are defined by (24) for free boundary conditions and by (pq) for periodic 
boundary conditions. Likewise the wavenumbers fa, fa, fa are given by (18. a) and by 
(18. b) for free- and for periodic boundary conditions, respectively. 
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